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INTRODUCTION 

Knowledge  of  the  neutron  intensity  in  the  core  is  needed  in 
many  applications  connected  with  the  operation  and  utilization  of 
a  nuclear  reactor.  This  intensity  is  defined  in  tejnns  of  the 
neutron  fltix  and  expressed  in  units  of  neutrons  per  squared  centi- 
meters per  second.   Much  has  been  written  describing  the  design 
and  the  capabilities  of  the  TRIGA  Mark  II  reactor  (26,  29,  30), 
but  very  few  results  have  been  published  describing  in  detail  the 
experimentally  measured  neutron  flux  distribution  in  the  TRIGA 
core.  G.  B,  West  (31)  of  General  Atomics  has  performed  multi- 
group  calculations  of  the  neutron  flux  in  the  homogenized, 
uniform  core.  These  calculations,  however,  frequently  do  not 
provide  adequate  information,  hence  it  is  desirable  to  measure 
experimentally  the  flux  in  the  reactor  core.   Experimental  veri- 
fication of  these  calculated  values  are  particularly  needed  for 
irradiation  facilities  such  as  the  central  thimble  and  the  rotary 
specimen  rack. 

Measurement  of  the  neutron  flux  distribution  and  calculation 
of  the  reactor  power  level  by  the  integrated  flux  method  in  the 
Kansas  State  TRIGA  core  was  suggested  by  W.  R.  Kimel,  head  of  the 
Nuclear  Engineering  Department  at  Kansas  State.  Work  on  this  pro- 
ject was  first  perfoinned  by  Robert  Ihde  using  indium  foils  as 
neutron  activation  detectors.   The  results  of  his  work  are 
available  in  the  form  of  an  incompleted  Master's  thesis  entit- 
led Neutron  Flux  Measurements  and  Power  Calibration  in  the 
Kansas  State  University  TRIGA  Mark  II  Nuclear  Reactor.  Special 


credit  must  be  given  to  Mr.  Ihde  for  his  efforts  in  procuring  the 
Reuter-Stokes  miniature  fission  chamber  used  for  the  flux  sceui- 
ning  in  my  work,  and  for  his  many  helpful  suggestions. 

This  work  attempts  to  provide  for  the  user  of  the  Kansas 
State  University  TRIGA  Mark  II  reactor  (17)  experimental  data 
about  the  neutron  flux  useful  for  planning  and  performing  experi- 
ments in  the  TRIGA  core.  The  specific  objective  of  the  work  is 
twofold:   first,  to  determine  experimentally  the  neutron  distri- 
bution in  the  Kansas  State  TRIGA  Mark  II  reactor  core  and  asso- 
ciated irradiation  facilities;  second,  to  estimate  the  reactor 
power  level  using  the  measured  distribution  in  the  core.  For 
neutron  activation  experiments  it  is  convenient  to  express  neu- 
tron flux  by  "Westcott  convention"  (32)  (defined  as  the  neutron 
density  times  the  velocity  of  2200  meters  per  second) .   Neutron 
flux  defined  in  this  manner,  together  with  the  cOjsorption  cross 
section  evaluated  at  2200  meters  per  second,  gives  for  a  1/v 
cibsorber  the  correct  absorption  rate  regardless  of  the  energy 
distribution  of  neutrons.  Occasionally,  for  activation  of 
resonance  detectors  or  fast  neutron  threshold  detectors,  addi- 
tional information  is  needed  about  the  neutron  velocity  distribu- 
tion in  the  core.  Values  are  therefore  est2±»lished  for  thermal, 
resonance,  and  fast  fluxes  in  two  of  the  most  frequently  used 
irradiation  facilities,  the  central  thimble  and  the  rotary  speci- 
men rack . 

For  measurement  of  the  absolute  neutron  flvix  in  the  TRIGA, 
standard  neutron  activation  techniques  (18,  22}  were  used.  The 
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determination  of  an  absolute  flux  is  in  general  difficult;  there- 
fore, it  was  decided  to  measure  the  fliix  carefully  at  a  single 
point  and  use  this  point  as  the  reference  for  all  other  flxix  meas- 
urements. For  convenience  in  defining  flxix,  the  point  should  be 
in  a  region  of  well  thermalized  flux.  Such  a  point  is  available 
in  the  thermal  column  of  the  TRIGA  reactor.   It  is  hoped  that 
calibration  of  flux  in  the  thermal  column  will  be  useful  not  only 
in  this  work,  but  as  a  standard  flux  value  for  future  experiments 
at  Kansas  State.  Two  independent  methods  were  employed  to  obtain 
the  eibsolute  activity  of  foils  irradiated  at  the  calibration 
point  in  the  thermal  column.  The  first  method,  beta-gamma  coin- 
cidence, is  well  known  and  has  been  reviewed  by  several  authors 
(3,  12,  24).   The  second,  the  "svmi  peak"  method,  is  a  relatively 
new  and  promising  technique  developed  at  the  Institute  for  Nuclear 
Physics  Research  in  the  Netherlands  (5,  6,  7,  8). 

Neutron  activation  methods  developed  in  the  literature 
(13,  22,  27,  28)  are  also  used  for  measuring  thermal,  resonance, 
and  fast  flux  values  for  the  central  thimble  and  the  rotary  speci- 
men rack.  Since  neutron  activation  is  generally  a  complex  and 
time  consimiing  method,  it  is  difficult  to  get  results  for  a  large 
number  of  points.  Thus,  for  the  purpose  of  a  detailed  scemning 
of  the  neutron  flux  in  the  core  a  simple,  compact  method  is 
desirable.   The  method  selected  utilizes  a  miniatiire  fission 
chamber  capa±)le  of  operating  over  large  ranges  of  neutron  flux 
with  a  high  degree  of  spatial  resolution  (1,  4,  25). 

The  total  power  generated  in  the  core  of  the  TRIGA  may  be 
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calculated  from  the  flux  plots  obtained  with  the  fission  chamber. 
This  approach  for  calculating  the  power  output,  called  the  "inte- 
grated flux"  method  (21,  22),  is  simple  in  principle;  but  it  is 
complicated  considerably  by  the  perturbations  introduced  as  a 
result  of  the  partially  inserted  control  rods  and  by  the  hetero- 
geneous nature  of  the  TRIGA  core. 

Section  (1.0)  deals  with  neutron  flxxx  measurements  in  the 
reactor  core  and  associated  irradiation  facilities.  It  is  divid- 
ed into  three  parts.  The  first  is  the  absolute  thermal  flux  meas- 
urement at  the  selected  point  in  the  thermal  colviinn.   The  second 
is  the  measurement  of  thermal,  resonance,  and  fast  fluxes  in  the 
central  thimble  and  the  rotary  specimen  rack.  The  third  is  the 
flux  scanning  with  the  miniature  fission  chamber.  Section 
(2.0)  presents  the  reactor  power  level  calculations  using  the 
flux  measurements  taken  in  Section  (1.0).  The  concluding  section 
summarizes  the  results  of  Sections  (1.0)  and  (2.0). 


(1.0)  Measurement  of  Neutron  Plux  in  the  Reactor  Core  and 
Irradiation  Facilities 

For  all  irradiations  2md  flux  measurement  in  the  Kansas 
State  TRIGA  Mark  II  core  and  irradiation  facilities,  Power  Level 
Channel  No.  1  on  the  reactor  console  was  used  as  the  reactor 
power  level  indicator.   This  channel  consists  of  a  compensated 
ion  chamber  that  feeds  a  raulti-range  linear  recorder  coupled  with 
an  automatic  controlling  device  for  maintaining  the  selected 
power  level  within  ±  1%  for  any  one  experiment  (17) .   Reactor 
power  levels  from  day  to  day  are  not  this  reproducible.  Experi- 
ence with  the  miniature  fission  chamber  and  a  BF3  ion  chamber 
introduced  in  the  central  thimble  indicates  that  standard  devia- 
tion between  identical  power  levels  over  relatively  short  periods 
(i.e.,  about  two  weeks)  is  from  ±  4%  to  ±  5%.   This  corapares  to 
a  value  of  ±  3%  quoted  by  General  Atomics  in  the  TRIGA  Operations 
and  Maintenance  Manual  (29) . 

Throughout  this  work  "Westcott  convention"  (32)  has  been 
used  to  define  neutron  flux.   In  this  definition  flux  is  written 
as  ♦  "  nv  where  the  quantity  n  is  the  total  neutron  density 
measured  by  the  neutron  detectors,  and  the  velocity  v^  is  the 
most  probable  speed  in  a  Maxwellian  distribution  at  a  temperature 
of  293*'k  (2200  meters  per  second) ,   This  method  of  defining  neu- 
tron flux  is  particularly  useful  in  experiments  requiring  the  cal- 
culation of  reaction  rates  because  it  eliminates  the  necessity  of 
knowing  the  mean  velocity  of  neutrons  in  the  core  and  because  for 
a  1/v  absorber  it  allows  the  use  of  the  cross  section  at  2200 
meters  per  second. 


6 

(1.1)  Measurement  of  Absolute  Neutron  Flux  in  the  Thermal  Colximn 

It  would  be  difficult  to  measure  neutron  flux  in  absolute 
units  at  every  position  required  in  these  experiments.  An  alter- 
native is  to  measxire  the  absolute  flux  at  a  selected  calibration 
point  and  use  this  standardized  flux  as  the  reference  for  other 
flxix  measurements.   It  is  important  to  choose  the  calibration 
point  at  a  position  where  the  velocity  distribution  of  neutrons 
is  known.   Since  "Westcott  convention"  is  used  to  define  flux, 
the  calibration  point  should  be  in  a  region  of  well  theirmalized 
flux  so  that  the  contribution  to  the  induced  activity  from  res- 
onance neutrons  is  small.   A  point  which  is  accessible  and  which 
satisfies  the  above  requirements  is  provided  in  the  TRIGA  ther- 
mal column.   (See  Fig.  1) 

All  irradiations  in  the  thermal  column  were  performed  at 
an  indicated  power  level  of  10  kilowatts  in  the  one  inch  by  one- 
quarter  inch  irradiation  groove  located  closest  to  the  center  of 
the  column.   The  groove  is  located  two  inches  above  the  horizon- 
tal centerline,  four  inches  right  of  the  vertical  center line,  and 
extends  three  feet  eight  inches  into  the  col\jmn.   The  foils  were 
positioned  four  inches  from  the  inner  end  of  the  groove  with  a 
wooden  stringer  normally  used  as  a  foil  holder  in  the  KSU  sub- 
critical  graphite  pile.   For  detailed  drawings  of  the  TRIGA  ther- 
mal colxomn  refer  to  the  TRIGA  Operation  and  Maintenance  Manual 
(29). 

The  irradiation  procedure  in  the  thermal  column  is  further 
complicated  by  high  radiation  levels  in  the  reactor  bay  if  the 
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thermal  col\ainn  door  is  open  at  a  power  of  10  kilowatts  (20) , 
Therefore  for  each  irradiation  the  sample  was  inserted  and  the 
thermal  column  door  closed  before  reactor  start-up.  The  reactor 
was  then  brought  to  a  power  of  10  kilowatts  on  a  period  of  less 
than  15  seconds.   Timing  began  the  instant  the  10  kilowatt  level 
was  reached.   The  maximum  additional  activity  introduced  by  this 
procedure  for  a  15  minute  irradiation  is  approximately  0.05%  of 
the  total  activity. 

By  measuring  the  absolute  initial  activity  of  a  suitable 
neutron  activation  detector  irradiated  in  the  thermal  column, 
neutron  flux  can  be  calculated  from  the  expression  (16) 

A«  =  VZja-e~^^)  ,  (1) 

o     a 

where 

V  =  the  volume  of  the  irradiated  detector, 

Z  =  the  macroscopic  activation  cross  section 
*   of  the  detector  isotope, 

^  =  the  average  neutron  flux  in  the  detector, 

X  =  the  radioactive  decay  constant  of  the 
detector  isotope,  and 

t  =  the  irradiation  time. 

Absolute  activities  are  determined  at  the  selected  position  using 

two  techniques  described  in  Sections  (1.1-a)  and  (1.1-b), 

(1.1-a)   Beta-Gamma  Coincidence  Method  for  Measuring  the  Activity 
of  Gold 

For  the  beta-gamma  coincidence  technique  a  neutron  activa- 
tion detector  that  emits  a  beta-ray  in  coincidence  with  a  gamma- 


ray  is  irradiated  in  the  thermal  column.   After  irradiation  this 
isotope,  decaying  at  the  rate  of  A  disentegrations  per  second, 
is  counted  with  two  detectors,  one  that  is  beta  sensitive  and 
one  that  is  gamma  sensitive.   The  count  rates  obtained  are 


^3  =  ^^B 


(2) 


where 


n  =  the  count  rate  of  the  beta  sensitive 
^   detector  and 

e  =  the  efficiency  of  the  beta  detector  for 
^   beta  radiation. 


and 


where 


n  =  Ae 
Y     Y 


(3) 


n  =  the  count  rate  of  the  gamma  sensitive 
^   detector  and 

e   =  the  efficiency  of  the  gamma  detector 
^   for  gamma  radiation. 

The  chance  of  observing  a  count  rate  in  both  detectors  simul- 
taneously is  given  by  e.e  and  the  true  coincidence  count  rate 
n.  is  given  by 


n.  =  Ae^E  , 
Py     3  y 


(4) 


By  combining  equations  (2),  (3),  and  (4),  the  following  expres- 
sion for  absolute  activity  is  obtained: 


n.n 
A  =  -^ 
^3Y 


(5) 
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In  practice  the  ideal  conditions  above  are  difficult  to 
satisfy?  thus  to  get  good  results  from  the  beta-gairana  coincidence 
method,  it  is  necessary  to  apply  various  correction  factors  (9, 
24)  to  equation  (5) .  For  the  particular  coincidence  arrangement 
used  here  the  only  corrections  considered  were  for  random  (acci- 
dental, chcince)  coincidences  and  for  background  coincidences. 
If  n*  is  the  total  observed  coincidence  rate,  if  n  is  the  ran- 
dom  coincidence  count  rate  due  to  finite  resolving  time  of  the 
coincidence  circuit,  and  if  n^^  is  the  background  coincidence  rate 
resulting  from  cosmic  coincidences,  then  the  true  covint  rate  can 
be  written  as  n^  »  n*  -  n^.  -  nj^. 

For  two  completely  independent  count  rates  (nj  and  n2) 
obtained  from  low  efficiency  detectors,  the  random  count  rate  is 
given  by  n  =  2Tnin2  where  t  is  the  resolving  time  of  the  coinci- 
dence circuit. 

A  gold  foil  one  mil  thick  and  one-half  inch  in  diameter  was 
used  as  the  neutron  detector  for  the  coincidence  determination. 
Gold  was  selected  for  two  reasons:   first,  the  decay  scheme  of 
gold  (24)  is  such  that  the  beta-gamma  method  applies;  second,  the 
activation  cross  section  and  other  necessary  material  constants 
are  well  established. 

The  block  diagram  in  Fig.  2  gives  the  particular  arrange- 
ment used  in  the  coincidence  counting  of  gold  performed  for  this 
work.  The  gamma  detector  shown  is  a  three  inch  by  three  inch 
Nal(Tl)  crystal  and  the  beta  detector  is  a  one  inch  by  one-quart- 
er inch  plastic  scintillation  crystal.  The  Tektronix  oscillo- 
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Figure  2.   Block  diagram  of  beta-gamma  coincidence  system, 


Equipment:   HV-1  John  Fluke  Power  Supply  Model  A0080A,  NE-179,  1300  volts. 
HV-2   Fluke  HV  Power  Supply  Model  412B,  NE-1198,  900  volts. 
Tennelec  TC  200  Amplifier. 

Tektronix  Type  545  Oscilloscope,  NE-28A .  « 

Tektronix  Type  53/54L  Plug-in-Unit,  NE-0284. 
Baird  Atomic  Scaler,  NE-145 . 

Beta  Detector  -  l"xl/4"  Plastic  Scintillation  Crystal,  NE-100. 
Gamma  Detector  -  3"x3"  Nal(Tl)  Crystal. 
TMC  4096  Multiparameter  System. 
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scope  in  the  coincidence  system  was  used  as  a  beta  pulse  dis- 
criminator and  shaper.   The  Tektronix  plug-in-unit  merely  inverts 
the  pulses  so  they  can  be  accepted  by  the  scaler.   The  coinci- 
dence gating  circuit  provided  in  the  TMC  4096  multiparameter 
system  may  be  operated  in  either  the  coincidence  or  anticoinci- 
dence mode.   The  gating  pulse  provided  by  the  beta  detector  must 
have  a  positive  voltage  between  three  and  twelve  volts,  and  it 
must  precede  the  gamma  pulse  by  250  nanoseconds.   Amplification 
and  shaping  delays  the  beta  pulse  making  it  necessary  to  intro- 
duce a  delay  line  in  the  gamma  branch  of  the  circuit.   Proper 
adjustment  of  the  delay  allows  the  beta  pulses  to  precede  the 
gamma  pulses  by  the  specified  time  interval. 

For  each  activity  determination  two  counting  periods  are 
necessary,  one  with  the  gating  circuit  in  the  coincidence  mode 
of  operation  and  the  other  with  the  gating  circuit  in  the  anti- 
coincidence mode.   Checking  the  coincidence  and  anticoincidence 
count  rates  for  each  deteinnination  affords  greater  confidence  in 
the  correct  performance  of  the  equipment. 

By  disregarding  background  and  random  coincidences ,  the 
fractional  standard  deviation  on  A  as  derived  by  Campion  and 
Taylor  (10)  is 

h 


^A  ^ 


-      (2e-e  -e„-e  +1) 
Hg^'   6  Y   6  y 


(6) 


Data  taken  for  the  gold  foil  irradiated  for  15  minutes  in 
the  thermal  colximn  are  summarized  in  Table  I.   Background  has 
been  sxibtracted  from  all  count  rates. 
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Table  I 
Beta-Gamma  Coincidence  of  Gold 


Time  After  n  n           A  o 

Removal       y                       B  BY  anticoinc.  A 

(hovirs)  (c/5  min.)  (c/5  min.)  (c/5  min.)  (c/5  min.)  (dps)  (dps) 

261.26  201,965     54,695  11,647  190,391    3,161  ±57 

355.52      71,942     22,463  4,744       67,160    1,135  ±32 

409.85      41,298     13,945  2,723      37,670      701  ±26 


For  the  above  case  the  number  of  random  coincidences  is  negligible 
compared  to  the  statistical  deviation  in  the  coincidence  count 
rate.   Initial  activity  A  is  calculated  from  the  activity  data 
using  the  weighted  least  squares  extrapolation  program  (Appendix 
A-1)  with  corrections  for  using  the  standard  deviations  from 
Table  I.   The  value  of  initial  activity  A^  is  corrected  for  flux 
depression  in  the  vicinity  of  the  foil  using  the  expressions 
given  in  the  Handbuch  der  Physik  (2) .   The  flux  depression  factor 
calculated  for  a  one  mil,  one-half  inch  gold  foil  in  a  graphite 
medium  is  1.0032.   The  cadmium  ratio  at  the  calibration  point  is 
79.5,  making  the  contribution  to  the  activity  from  the  resonance 
flux  region  negligibly  small. 

The  corrected  value  for  initial  saturated  activity  of  the 
gold  foil  is  (2.60  ±  0.07)  x  10**  dps/gm-watt.   By  using  this 
initial  activity  and  98.7  ±  0.6  barns  as  the  activation  cross 
section  of  gold,  the  thermal  neutron  flux  at  the  selected  cali- 
bration point  in  the  thermal  column  is  (8.80  ±  0.25)  x  10 •♦ 
n/cm^- sec-watt. 
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(1.1-b)   The  "Sum  Peak"  Method  for  Determining  the  Activity 
of  Scandixam 

In  the  past  few  years  a  group  from  the  Institute  for 
Nuclear  Physics  Research  in  Amsterdam,  Holland  has  developed  a 
means  of  determining  the  absolute  activity  of  a  radioactive 
source  which  they  call  the  "s\m   peak"  method.   In  this  method, 
an  isotope  emitting  two  gamma- rays  in  coincidence,  each  100% 
abundant,  is  counted  using  a  scintillation  crystal  and  a  multi- 
channel analyzer.  A  typical  gamma  energy  spectrum  for  this  case 
is  shown  in  Fig.  3. 

In  this  figure,  Ai  and  A2  denote  the  areas  under  the  two 
photopeaks,  A12  is  the  area  \inder  the  sum  peak,  and  T  is  the 
total  area  under  the  spectrum.   Now,  defining  the  photoelec- 
tric efficiencies  ei  and  e2  as  the  probabilities  that  a  gamma-ray 
is  captured  in  the  crystal  by  a  photoelectric  process,  the  total 
efficiencies  ti  and  t2  as  the  probabilities  that  a  gamma-ray  is 
captured  in  the  crystal  by  any  means,  and  the  activity  A  as  the 
total  number  of  disintegrations  in  the  sample,  the  following  equa- 
tions can  be  derived  (5) : 

Aj  =  Aei (l-t2)  /  (7) 

A2  =  Ae2(l-ti) ,  (8) 

A12  =  Aeie2f  (9) 

T  =  A(ti  +  t2  -  titz) .  (10) 

Combining  the  four  equations  above  yields: 

A  =  A1A2/A12  +  T.  (11) 

Thus  with  a  single  measurement  of  the  gamma  spectriim  the  abso- 
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Figure  3.   Gamma  spectrum  of  Sc 
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lute  activity  can  be  determined. 

Scandium  was  chosen  as  the  activation  detector  for  the  "sum" 
peak"  technique  because  it  exhibits  only  two  cascading  gammas. 
Methods  for  treating  more  complex  decay  schemes  by  the  "sum  peak" 
method  are  being  developed,  but  to  simplify  the  calculations 
it  was  convenient  to  use  scandium  with  its  simple  decay  scheme. 

Estimates  of  the  errors  involved  in  the  three  activity 
determinations  by  the  "simi  peak"  method  indicate  an  uncertainty 
of  approximately  ±  13%  for  each  determination.   Two  difficulties 
in  the  method  result  in  this  relatively  high  error.   The  first 
is  the  uncertainty  in  calculating  areas  Aj ,  A2,  and  A12  so  that 
they  correspond  to  the  quantities  represented  by  equations  (7) , 
(8) ,  and  (9) .   Area  Aj  is  particularly  difficult  because  the 
Compton  contribution  from  gamma  #2  is  superimposed  on  the  spec- 
trum in  this  energy  region.   The  second  difficulty  results 
because  in  practice  the  first  few  channels  of  the  analyzer 
cannot  be  used.   This  means  that  to  integrate  the  entire  spectriim, 
an  extrapolation  to  zero  energy  in  required.   This  extrapolation 
may  introduce  appreciable  inaccuracies  in  T, 

Scandium  in  the  form  of  powdered  SC2O3  was  irradiated  for 
30  minutes  at  the  calibration  point  in  the  thermal  column,  and  the 
gamma  energy  spectrum  was  obtained  using  the  TMC  4096  Multipara- 
meter system.   Numerical  integrations  were  performed  to  de term- 
mine  areas  Ai ,  A2,  Ai2»  and  T  from  this  scandium  spectrum.   Data 
obtained  from  these  integrations  are  given  in  Table  II.   Flux 
values  have  been  corrected  for  flux  depression  of  the  sccindium 


17 
sample  in  graphite.  The  activation  cross  section  for  scandium 
used  in  the  flux  calculation  is  24  ±  0.1  barns  (19),  and  the 
average  value  of  neutron  flux  obtained  by  the  "sum  peak"  method 
was  (10.47  ±  1.36)  x  lO**  n/cm^- sec-watt. 
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(1.2)   Measurement  of  Thermal,  Resonance,  and  Fast  Fluxes  in  the 
Central  Thimble  and  the  Rotary  Specimen  Rack 

Since  the  central  thimble  and  the  rotary  specimen  rack  are 
commonly  used  for  a  variety  of  experiments  on  the  TRIGA,  it  is 
desirable  to  provide  for  the  TRIGA  user,  detailed  information 
about  the  neutron  flux  in  these  facilities.   In  a  thermal  reactor 
neutrons  are  continuously  born  in  the  fission  process  with  ener- 
gies from  0.05  to  15  iMev.   The  fission  neutrons  are  slowed  down 
in  the  moderator  material  until  they  reach  a  state  of  thermal 
equilibrium  with  the  moderator  molecules.   Figure  4  shows  the 
approximate  shape  of  the  neutron  energy  spectrxim  for  a  thermal 
reactor  and  gives  the  boxindaries  of  the  various  flux  regions. 

Gold,  indium,  and  sulfur  induced- activation  detectors  were 
irradiated  in  the  central  thimble  and  the  rotary  specimen  rack 
using  the  standard  polyethylene  sample  holders  provided  for  this 
purpose  (29).   The  detectors  were  positioned  1/8  inch  from  the 
bottom  of  these  holders  and  two  1/8  inch  aluminxim  disks  designed 
to  hold  the  detectors  in  place.  For  irradiations  in  the  central 
thimble  the  holders  were  filled  with  water  to  prevent  perturba- 
tions in  the  flux  due  to  an  air  gap  not  ordinarily  present  in  the 
core  at  this  position.   Detectors  irradiated  in  ttie  central 
thimble  were  positioned  14  inches  below  the  upper  edge  of  the  top 
grid  plate.   This  position  represents  the  point  of  maximum  flux 
in  the  core  and  is  the  normal  position  for  irradiations  performed 
in  the  central  thimble.   For  irradiations  in  the  rotary  specimen 
rack  the  sample  holders  were  simply  lowered  into  the  position 
directly  below  the  insertion  tube. 
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Figure  4.   Representative  neutron  energy  spectrum  for 
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Figure  5.   Block  diagram  of  the  Geiger-Mueller  counting  system. 
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Accurate  removal  and  insertion  times  were  obtained  using  a 
Standard  Electric  Timer  (Type  SM-60,  NE-88) .   Upon  insertion  the 
time  was  noted,  and  after  the  desired  amount  of  irradiation  time 
the  sample  was  pulled  from  the  flux  region.   The  maximum  esti- 
mated error  in  this  procedure  is  ±  1.0  seconds  for  both  removal 
and  insertion;  therefore,  the  maximum  error  in  irradiation  time 
is  approximately  2.0  seconds. 

After  removal  from  the  reactor  the  samples  were  counted 
with  the  Geiger-Mueller  counting  system  represented  by  the  block 
diagram  in  Fig.  5.   Counting  data  obtained  with  this  system  were 
used  to  calculate  the  initial  saturated  count  rate  for  each  de- 
tector.  This  is  accomplished  by  using  the  weighted  least  squares 
extrapolation  program  described  in  Appendix  A-1.   The  Geiger- 
Mueller  system  is  calibrated  for  measurement  of  absolute  flux  by 
counting  gold  foils  irradiated  in  the  standard  flux  in  the  ther- 
mal  column. 

(2.2-a)   Thermal  Flux 

Gold  and  indium  were  chosen  as  activation  detectors  for 
measurement  of  thermal  flux  because  their  activation  cross  sec- 
tions are  well  known,  and  because  1/v  energy  dependence  can  be 
assumed  for  both  materials  in  the  thermal  region.   Both  gold  and 
indium  exhibit  large  resonance  absorption  peaks  in  the  slowing 
down  region>  thus  to  separate  thermal  activity  from  resonance 
activity,  bare  and  cadmium  covered  foils  were  irradiated  in  each 
position.   Cadmium  20  mils  thick  has  a  high  capture  cross  section 
for  thermal  neutrons;  however,  the  cross  section  drops  abruptly 
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for  neutrons  with  energies  greater  than  0.5  ev.   Consequently  the 
bare  foil  measures  both  thermal  response  and  resonance  response; 
whereas,  the  cadmium  covered  foil  measures  only  resonance  neutrons, 

At  this  point  two  correction  factors  need  to  be  considered. 
One,  Fpj»  accounts  for  the  epithermal  absorption  of  neutrons  in 
cadmium;  the  other,  Fj.u»  corrects  for  flux  depression  in  the 
vincinity  of  the  foils.   Procedure  for  determining  these  correc- 
tion factors  is  treated  in  detail  in  the  Handbuch  der  Physik  (2) . 
Values  for  these  corrections  applied  to  detectors  used  in  this 
work  are  summarized  in  Table  III. 

Table  III 


xs=c:srac:s 


Correction  Factors  For  Activation  Detectors 


Detector 


Thickness  Radius  Surrounding 


(mil) 


(inches)   Medium 


cd 


th 


Gold 

1 

1/2 

Water 

1.00 

1.0168 

Gold 

1 

1/2 

Graphite 

1.00 

1.0032 

Gold 

5 

1/2 

Air 

1.00 

1.0000 

Gold 

5 

1/2 

Graohite 

1.00 

1.0150 

Indivm 

2 

1/2 

Water 

1.05 

1.0342 

Indium 

2 

1/2 

Air 

1.00 

1.0000 

Indiian 

2 

1/2 

Graphite 

1.05 

1.0035 

Scandixim 

10 

1/2 

Graphite 

1.0035 

Applying  the  above  correction  factors  and  subtracting  the 
activity  of  the  cadmium  covered  foil  from  the  activity  of  the  bare 
foil  yields  the  following  equation  for  the  activity  due  to 
thermal  neutron  flux  (2) : 


A^.,  =  F. ,  (A,  -  F^  ,A^  ,)  , 


'th 


th'"b 


cd  cd' 


(12) 
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where 

A  .  ■  the  activity  in  the  foil  due  to  thermal 
neutrons , 

A-   -  the  activity  of  the  bare  foil,  and 

A  .  "  the  activity  of  the  cadmium  covered  foil. 

(2.2-b)   Resonance  Flux 

The  same  gold  and  indium  foils  used  for  measurement  of 
thermal  flux  were  also  used  for  the  measurement  of  the  slowing 
down  density  at  specific  energies  in  the  resonance  region.  The 
single,  large  resonance  peak   exhibited  by  both  detectors  accounts 
for  roost  of  the  activation  in  the  epi thermal  region,  giving  an 
indication  of  the  neutron  flux  at  the  resonemce  energies.   The 
cadmium  difference  method  (22)  was  then  used  to  establish  the 
1/E  dependence  of  the  spectrum.   The  cadmium  ratio  is  expressed 
as 

R  »  1  +    th  (13J 

^cd 

For  a  1/E  energy  spectrum  the  ratio  ^j./^^   will  be  approximately 
constant.   Ruzik  (22)  derives  the  following  expression  for  this 
ratio: 

♦r/*th  "  (R-1) (1+a)  '  ^^^^ 

where  a  is  a  constant  for  each  material  given  by 

/o^dE/E 


.  (15) 

/o^/^dE/E 
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If  the  spectrum  is  truly  1/E,  the  differential  flux  is  then  given 
as 

♦(E)  =  ♦^/E  .  (16) 

(2.2-c)   Fast  Flux 

Sulfur  was  chosen  as  the  threshold  detector  for  measurement 
of  the  integrated  fast  flux  in  the  high-energy  part  of  the  neutron 
spectrum.   The  reaction  of  interest,  S^2(n,p)p32,  has  a  saturation 
cross  section  of  300  millibarns  and  an  effective  threshold  energy 
of  2.9  Mev.  (11).   The  threshold  energy  of  sulfur  is  established 
by  representing  the  sulfur  cross  section  as  an  idealized  step 
function.   This  step  function  is  zero  for  all  neutron  energies 
below  2.9  Mev.  and  300  millibarns  for  energies  above  2.9  Mev. 
(See  Fig.  6) 

The  sulfur  detectors  were  pressed  into  pellets  3/16  inch 
thick  and  3/4  inch  in  diameter.   Approximately  three  grams  of  USP 
grade,  powdered  sulfur  was  pressed  with  no  binder  and  no  heating 
at  an  indicated  pressure  of  2000  psi  on  a  hydraulic  press.  A 
phosphorous  pellet  of  the  same  dimensions,  3/16  inch  thick  and 
3/4  inch  in  diameter,  was  needed  for  determining  the  efficiency 
of  the  Geiger-Mueller  counting  system  for  the  P^^  betas.  Amor- 
phorous  powdered  red  phosphorous  was  pressed  at  a  temperature  of 
400  C  and  a  pressure  of  2000  psi,  using  approximately  four  weight 
percent  polystyrene  binder.   After  pressing,  the  pellets  were 
sprayed  with  Krylon  to  prevent  flaking. 

For  calibration  of  the  Geiger-Mueller  system  the  phospho- 
rous pellet  was  irradiated  in  the  known  thermal  flux  of  the  ther- 
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Figure  6.   Fission  spectrum  and  cross  section  for 
S^2(n,p)P'''  reaction. 


Effective  threshold  energy  =  2.9  Mev 
Activity  of  the  sulfur  detector  =  Area  A  +  Area  C 
Integrated  fast  flux  above  2.9  Kev  =  Area  B  +  Area  C 
Area  A  =  Area  B 
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mal  column  and  counted  with  the  Geiger-Mueller  system.   Since  the 
remge  of  the  P^^  betas  (E.  =  1.70  Mev.)  is  considerably  less  than 
the  thickness  of  the  pellets  and  since  the  density  of  the  two 
materials  is  nearly  the  same,  it  is  assumed  that  the  active  weight 
seen  by  the  counting  system  is  the  same  for  phosphorous  as  for 
sulfur.   With  this  assumption  the  following  expression  can  be 
obtained  for  fast  flux  (22) : 

Vs^s^p'^th  .  (17) 


where 


*f  "   K  C  A  O 
s  p  p  o 


()>,  =  the  integrated  fast  flux  above  2.9  Mev., 

S  =  the  content  of  P^^  in  the  phosphorous 
P   pellet  in  weight  percent, 

K  =  the  isotopic  content  of  S'^, 

s 

C  =  the  initial  saturated  count  rate  from  the 
P   phosphorous  pellet, 

C  =  the  initial  saturated  count  rate  from  the 
^   sulfur  pellet, 

A  =  the  atomic  weight  of  sulfur, 
s 

A  =  the  atomic  weight  of  phosphorous, 

o  =  the  saturation  cross  section  for  the 
°   S^2(n,p)P^^  reaction, 

o  =  the  activation  cross  section  for  P^^,  and 
P 

«|.  .  =  the  theiroal  flux  at  the  calibration  point 
in  the  thermal  column. 


(2.2-d)   Tabulated  Results  for  the  Central  Thimble  and  the 
Rotary  Specimen  Rack 

Table  IV  lists  the  constants  of  each  detector  used  in  the 

flux  calculations  (14,  19,  22).   The  thermal  activation  cross 
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sections  are  given  for  a  neutron  velocity  of  2200  meters  per 
second.   Alpha  is  the  ratio  of  the  resonance  component  of  the 
resonance  activation  to  the  1/v  component  and  is  defined  by 
equation  (15) .   Table  V  siommarizes  irradiation  and  counting  data 
for  foils  used  in  the  flux  determinations.   The  detectors  irra- 
diated in  the  thermal  column  at  the  flux  calibration  point  were 
used  for  the  calibration  of  the  Geiger-Mueller  system  counting 
system.   The  results  of  this  calibration  for  gold  foils  are  given 
in  Table  VI.   Table  VII  shows  the  results  of  the  cadmixom  differ- 
ence method  for  gold  and  indium  detectors.   The  ratio  of  <^   Z+th' 
is  approximately  constant  in  the  central  thimble  and  the  rotary 
specimen  rack,  indicating  that  a  1/E  spectrum  is  present. 

Table  IV 


Detector  Constants 

Detector        Reaction  Activation  a 

Cross  Section 
(barn) 

Gold  Aul57   n  +  Aul37  ^  j^^lSQ  93.7  35.4 

(2.7  day)     3  =  .960  Mev.  (98.6%) 
Y  =  .412  Mev.  (100%) 

Indium       n  +  In^^s  ^  m^l^  155  37.0 

ml  15 

(54.2  min.) 

Sulfur       n  +  s32->.p32+p+   Saturated 

S^^  Reaction  Cross 

(14.3  day)    6  =  1.70  Mev.  (100%)  Section  =0.30 

Phosphorous  n  +  P^^  •*  P^^  0.23 

p31 

(14.3  day)    B  =  1.70  Mev.  (100%) 
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Table  V 

===:=:=======: 

Irradiation  and  Counting  Data 

Indicated 

Initial 

Reactor 

Saturated 

Power   Irradiation 

Counting 

Level      Time 

Weight  Thickness  Rate 

Detector 

Position 

(watts)    (minutes) 

(grams) 

(mil 

s)  (cpm/gm) 

Au 

C.  Th. 

10      20 

.0696 

1 

.2024  X  10^ 

Au* 

C.  Th. 

10      20 

.0695 

1 

.5830  X  10^ 

Au 

RSR 

10      20 

.3496 

5 

.8027  X  10^ 

'J 

Au* 

RSR 

10      20 

.3514 

5 

.1138  X  107 

In 

C.  Th. 

10       3 

.0135 

2 

.5591  X  109 

In* 

C.  Th. 

10       3 

.0185 

2 

.1465  X  10^ 

In 

RSR 

10       3 

.0154 

2 

.9237  X  10^ 

In* 

RSR 

10       3 

.0160 

2 

.1626  X  10^ 

S 

C.  Th. 

10,000       3 

.1047  X  10^ 

_     o 

S 

C.  Th. 

1,000      20 

.1012  X  10^ 

S 

RSR 

10,000      40 

.3694  X  10^ 

S 

RSR 

90,000       5 

.3116  X  10^ 

Au 

Th.  Col. 

10,000      15 

.0695 

1 

.4228  X  10^ 

Au 

Th.  Col. 

10,000      10 

.3514 

5 

.1107  X  109 

P 

Th.  Col. 

10,000      30 

.1344  X  10^ 

Table  VI 

Calibration  of  the  Geiger-Mueller 

Counting  System 

Inverse  Counter 

Method  of  Absolute 

Efficiency 

Detectoi 
Gold  {1 

mil)    e 

Flux  Determination 

(dps/cpm) 

-  Y  Coincidence  Method 

(Aul9  8) 

37.7 

Gold  (1 

mil)    "1 

Sum  Peak"  Method  (Sc**^) 

43.5 

Gold  (5 

mil)     e 

-  Y  Coincidence  Method 

(Aul98) 

144 

Gold  (5 

mil)    "1 

Sum  Peak"  Method  (Sc**^) 

166 

*    Indicates  cadmium  covered  foils. 

> 

Table  VII 


Cadmium  Difference  Method  for  the 
Central  Thimble  and  the  Rotary  Specimen  Rack 
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Position 

Detector 

(R) 

*r/*th 

C.  Th. 

Gold  (Au^^®) 

3.46 

.0264 

C.  Th. 

Indium  (In^^S) 

3.63 

.0238 

RSR 

Gold  (Aul9  8) 

7.06 

.0107 

RSR 

Indium  (In^^S) 

5.42 

.0141 
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(1.3)   Neutron  Flux  Mapping  in  the  Reactor  Core 

The  spatial  neutron  distribution  in  the  TRIGA  core  was  ob- 
tained by  using  a  miniature  fission  chamber  as  a  neutron  detector. 
Access  to  the  core  is  available  through  the  0.314  inch  foil  inser- 
tion holes  (29)  provided  in  the  top  grid  plate.   (See  Fig.  7)   A 
fission  chamber  capable  of  fitting  into  these  foil  insertion 
holes  and  of  measuring  flux  in  the  range  from  1  x  10^  to  2  x  lO^** 
n/cm^-second  is  available  commercially  from  Reuter-Stokes 
Electronics  Components,  Inc.  of  Cleveland,  Ohio.   The  Reuter- 
Stokes  chamber  has  a  sensitive  length  of  1.0  inch  and  an  outside 
diameter  of  0.250  ±  0.005  inches.   The  inner  electrodes  are 
titanium  eind  the  outer  shell  is  304  stainless  steel.   The  neutron- 
sensitive  portion  contains  a  total  of  3.34  milligrams  of  uranium, 
93%  enriched  in  U-235,  surrovinding  an  ion  chamber  filled  with 
helium  gas  at  a  pressure  of  76  cm.  Hg. 

Twenty  feet  of  coaxial  cable  connects  the  chamber  with  the 
associated  equipment  located  on  top  of  the  reactor  shield.   A 
seal  between  the  active  volume  of  the  chamber  and  the  cable  pre- 
vents gas  migration  between  the  detector  and  the  cable,  providing 
uniform  sensitivity  over  temperatures  ranging  as  high  as  200  C. 

Figure  8  shows  the  physical  arrangement  for  positioning 
the  fission  chamber  in  the  core.   A  26  11/16  inch  rigid,  thin- 
walled  aluminxom  probe  containing  the  fission  chamber  slides 
through  the  foil  insertion  holes  in  the  top  grid  plate.   The 
small  hole  in  the  bottom  of  the  probe  allows  the  entire  asseit±)ly 
to  fill  with  water.   A  flexible  20  foot  long  aluminum  guide  tube 
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Figure  7.   Core  loading  diagram 

showing  the  foil  inser- 
tion holes. 


r   J    Fuel  element 
|fTr[h  Graphite  element 
O   Foil  insertion  hole 
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Figure  8.   Positioning  of  the  fission  chamber  in 
the  core 
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connects  the  probe  to  the  traversing  mechanism  on  the  top  of  the 
reactor  shield.   The  flexibility  of  the  guide  tube  allows  access 
to  all  the  foil  insertion  holes  except  I  and  J  despite  various 
obstructions  in  the  pool.   Holes  I  and  J  are  blocked  by  the  ter- 
minus assembly  of  the  pnevimatic  transfer  system. 

For  each  vertical  traverse  the  probe  containing  the  fission 
chamber  was  completely  inserted  as  shown  in  Fig.  8.   Since  the 
position  of  the  chamber  within  the  probe  is  fixed  and  the  length 
of  the  probe  is  known,  the  position  of  the  chamber  with  respect 
to  the  upper  surface  of  the  top  grid  plate  is  established.  With 
the  bottommost  position  as  zero,  the  probe  was  pulled  upward 
through  the  core  at  the  rate  of  one  inch  per  minute  using  the 
traversing  mechcinism.   Current  readings  were  taken  with  a 
Victoreen  micro-micro  ampmeter  at  0.50  inch  intervals.  During 
each  traverse  a  BF3  ion  chamber  positioned  two  inches  above  the 
top  grid  plate  in  the  central  thimble  provided  an  additional  check 
on  the  stability  and  reproducibility  of  the  reactor  power  level. 

(1.3-a)   Calibration  of  the  Fission  Chamber  for  Measurement  of 
Absolute  Neutron  Flxix 

To  calibrate  the  fission  chamber  for  the  measurement  of 

absolute  fl\ix  the  chamber  response  is  determined  in  the  region 

of  known  thermal  flux  in  the  thermal  column.   Table  IX  gives  the 

results  of  this  calibration. 
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(1.3-b)   Discussion  of  the  Neutron  Flux  Plots  in  the  Core 

Figures  9  through  23  are  vertical  plots  of  the  flux  in  the 
central  thimble  and  the  14  flux  insertion  holes  (A  through  P, 
excluding  I  and  J) .  Data  for  the  plots  were  taken  at  a  reactor 
pov/er  level  of  100  watts.   This  power  level  is  indicated  by  a 
reading  of  0.33  x  3  x  10 ^  watts  on  the  linear  recorder  on  the 
reactor  console.   All  flux  values  measured  with  the  fission  cham- 
ber were  corrected  for  flux  perturbations  caused  by  insertion  of 
the  chamber  into  the  core. 

Figures  24  and  25  show  the  response  of  the  miniature  fission 
chamber  as  the  reactor  power  level,  indicated  by  the  linear  rec- 
order, is  varied  from  10  watts  to  90  kilowatts.   The  linear  rec- 
order (29)  has  a  linear  scale  from  0  to  11  and  a  range  selector 
which  determines  the  value  of  the  full  scale  reading  (full  scale 
here  refers  to  a  reading  of  10  on  the  linear  scale) .   Figure  24 
shows  the  fission  chamber  response  as  the  power  indicated  by  the 
linear  recorder  is  increased  over  the  full  scale  (0  to  10)  at  a 
single  range  setting.   Note  that  the  response  is  linear  with  a 
slope  of  0.95  and  an  intercept  of  about  0.03  of  the  full  scale 
reading.   The  data  shown  in  Fig.  24  were  taken  at  a  range  selector 
position  of  10  x  10^  watts  but  a  similar  response  was  noted  for 
each  of  the  other  range  selector  positions.   Thus  the  power  indi- 
cated by  the  fission  chamber  for  any  linear  recorder  position  can 
be  read  directly  from  Fig.  24.  For  example  if  the  linear  rec- 
order reads  300  watts  (i.e.  0.30  of  full  scale  with  a  range  sel- 
ector position  of  10  x  10^  watts) ,  the  power  indicated  by  the 
fission  chamber  from  Fig.  24  is  0.315  of  10  x  10^  watts  or  315 
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watts.  Figure  25  shows  the  fission  chamber  response  over  the 
entire  power  range  (10  watts  to  90  kilowatts) .   The  data  at 
10  watts,  100  watts,  1  kilowatt,  and  10  kilowatts  were  taken  at 
0.33  of  full  scale  with  the  range  selector  in  the  appropriate 
position.  The  data  at  40  and  90  kilowatts  were  taken  at  0.40 
and  0,90  of  full  scale  respectively.  For  the  data  at  0.33  of  full 
scale,  the  slope  of  the  response  curve  is  1.0.  The  fission  cham- 
ber response  at  40  and  90  kilowatts,  however,  indicates  a  power 
reading  that  falls  below  the  extrapolated  line  of  slope  1.0. 
The  slope  of  the  response  curve  decreases  at  40  euid  90  kilo- 
watts because  the  slope  of  the  fission  chamber  response  versus 
linear  recorder  power  is  less  than  1.0  as  the  position  of  the 
recorder  pen  on  the  linear  scale  is  increased  (See  Fig.  24) . 

Throughout  this  work,  fission  chamber  flux  measurements  were 
made  at  a  recorder  pen  setting  of  0.33  of  full  scale.   The  only 
exceptions  to  this  procedure  were  measurements  taken  at  a  linear 
recorder  power  level  of  90  kilowatts  (See  Fig.  28  and  29),  in 
which  case  the  pen  recorded  0.90  of  full  scale. 

At  low  power  levels  (i.e.  <  10  watts)  the  gamma  radiation 
resulting  from  the  build  up  of  decaying  fission  products  may 
contribute  significantly  to  the  response  of  the  fission  chamber. 
The  magnitude  of  this  effect  depends  on  the  operating  history  of 
the  reactor,  but  in  general  for  power  levels  of  100  watts  or 
greater  it  is  small  compared  to  the  detector  response  to  neutrons. 
The  edjove  effect  is  not  shown  in  Pig.  25,  but  for  similar  data 
taken  after  long,  high  power  operation  of  the  reactor,  the  fission 
chamber  indicated  a  significantly  higher  reading  at  a  linear 
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recorder  power  of  10  watts. 

The  shapes  of  flux  plots  are  not  changed  as  the  power  level 
is  increased,  except  for  holes  M  and  N.   The  shape  change  in 
these  two  holes  is  a  result  of  the  change  in  the  regulating  rod 
position  as  the  power  level  is  increased  above  10  kilowatts. 
Figures  26  and  27  are  representative  plots  of  the  flvix  to  power 
ratio  for  traverses  at  100  watts  and  at  90  kilowatts  in  hole  M 
and  in  the  central  thimble.   A  result  similar  to  the  plot  in  M 
can  be  obtained  for  hole  N,  and  results  similar  to  the  plot  in 
the  central  thimble  can  be  obtained  for  the  other  foil  insertion 
holes.   Figures  28  and  29  further  indicate  the  relation  of  the 
flux  shape  to  the  power  level  by  plotting  the  fractional  dif- 
ference between  the  fliix  to  power  ratio  at  100  watts  and  90  kilo- 
watts versus  vertical  position. 

The  flux  shapes  are  consistent  with  what  might  be  expected 
for  the  TRIGA  core.   The  entire  distribution  seems  to  be  shifted 
slightly  downward  as  a  result  of  the  partially  inserted  neutron 
cdssorbing  control  rods  in  the  upper  portion  of  the  core.  Sharp 
flux  depressions  are  noticed  locally  around  the  control  rods  and 
definite  peaking  is  noticed  in  the  water  columns  below  the  control 
rods  and  in  the  central  thimble.   Flux  peaking  is  also  apparent 
in  the  graphite  end  reflectors.   The  samarium  oxide,  burnable 
poison  discs  located  between  the  fuel  and  the  end  reflectors, 
tend  to  depress  the  flux  in  this  region,  causing  even  more  pro- 
nounced flux  peaking  in  the  end  reflectors .   It  should  be  noted 
that  the  flux  peaks  in  the  plots  obtained  with  the  fission  cham- 
ber may  be  appreciably  larger  and  sharper  if  corrections  for  the 
finite  resolution  of  the  chamber  are  considered. 
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(2.0)   Power  Calibration 

Measurement  of  the  quantity  of  heat  produced  per  unit  time 
in  the  core  is  essential  to  the  safe  control  and  operation  of 
the  reactor.   Several  power  calibrations  have  been  performed  on 
the  Kansas  State  TRIGA  using  calorimetric  methods  (i.e.,  measure- 
ment of  the  increase  in  water  temperature  in  the  reactor  pool) . 
These  calorimetric  methods  are  adequate;  however,  measurement  of 
the  neutron  flux  distribution  in  the  core  provides  an  additional 
means  for  calculating  reactor  power  level. 

Each  fission  of  a  U-235  atom  .releases  energy  equivalent  to 
197  ±  8  Mev.  or  (3.16  +  0.13)  x  lO"^^  watt-sec  with  the  energy 
distribution  given  in  Table  X  (21) . 

Table  X 


Fission  Energy  Distribution 

Kinetic  energy  of  fission  fragments 165  ^+  5  Mev 

Prompt  Y~J^sys 

Kinetic  energy  of  fission  neutrons   5 

Fission  product  decay  )-  -^ 

Gamma  . 6 

Beta 5 


Neutrino 


Total 


11 


197  ± 


This  does  not,  however,  represent  the  total  heat  energy  released 
in  the  reactor  core  region.   Some  of  the  fission  energy  is  not 
converted  to  heat  in  the  core;  instead  it  escapes  the  core  region 
in  the  form  of  radiation.   Further,  some  additional  heat  results 
from  capture  gammas  formed  in  the  core.   These  effects  depend  on 
the  reactor  geometry  and  design  and  are  extremely  difficult  to 
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calculate.  Therefore,  in  this  work  the  total  power  generated  in 
the  reactor  is  taken  as  the  total  nuclear  energy  released  in  the 
fission  process,  specifically 


P  =  KVZ^4i 


(18) 


where 


P  =  the  total  power  output  of  the  reactor  in 
watts , 

V  =  the  volume  of  U-235  in  the  core, 

Z  =  the  macroscopic  fission  cross  section  for 
U-235  at  2200  meter/second, 

^  =  the  average  flux  in  the  fuel  region  of  the 
core ,  and 

K  =  the  total  nuclear  energy  released  per 
fission  in  watt-sec/fission. 

This  value  of  total  nuclear  power  cannot  be  compared  directly  to 
the  calorimetric  power  generated;  hov/ever,  by  simply  adjusting 
the  value  of  K  an  estimate  of  calorimetric  power  output  can  be 
obtained.   The  value  of  K  used  to  estimate  the  calorimetric  power 
was  (2.88  +  0.20)  x  10~^^  watt-sec  per  fission. 

By  assuming  that  the  fission  cross  section  is  known  for  a 
given  quantity  of  U-235,  the  problem  is  reduced  to  calculating 
the  average  flux  in  the  fuel  region  of  the  reactor  core.   The  ex- 
pression for  the  average  flux  is 


■il. 


<^dV  (19) 

V 

c 
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where 

V  =  the  total  volijine  of  the  fuel  in  the 
^       reactor  core,  and 

^   =  the  neutron  flux  as  a  function  of 
position  in  the  reactor. 

The  above  integration  is  performed  numerically  by  dividing  the 

core  into  small  volume  elements  as  shovm  in  Fig.  30.  Each 

annular  volume  element  is  one  inch  thick  and  includes  one  fuel 

element.  The  expression  for  ^  then  becomes 

M    N 


*  =  HN  E  E*ij  ^''^ 


.  i=i  j=i 

where  th 

♦• .  =  the  average  flux  in  the  ij   element, 

M  =  the  number  of  vertical  levels,  and 

N  =  the  number  of  fuel  elements. 

(3.1)   Ratio  of  the  Maximum  Flux  in  the  Moderator  to  the  Average 
Flux  in  the  Fuel 

In  a  heterogeneous  reactor  such  as  the  TRIGA  Mark  II  it  is 
not  practical  to  measure  directly  the  neutron  flux  within  the 
fuel  element.   Instead  the  flux  measurements  were  made  with  the 
miniature  fission  chamber  in  the  water  region  between  the  fuel 
elements.   From  these  measurements  the  average  flux  in  a  fuel 
element  for  a  particular  horizontal  level  is  obtained  by  averag- 
ing the  flux  values  at  the  four  corners  of  an  annular  element  and 
by  multiplying  this  average  by  the  ratio  of  the  maximum  fl\ix  in 
the  water  moderator  to  the  average  flux  in  the  fuel.   This  ratio 
is  approximately  equal  to  the  disadvantage  factor  c  and  can  be 
estimated  as  the  absorption  cross  section  of  the  hydrogen  in  the 
water  divided  by  the  absorption  cross  section  of  the  hydrogen  in 


Figure  30.   Loading  diagram  showing  the 
annular  elements.   (Note  that 
the  core  was  divided  similarly 
for  each  of  the  14  different 
vertical  levels.) 


C    j       Fuel  Element 
(f||]|)   Graphite  Element 
Q    Foil  insertion  hole 
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the  fuel.  By  using  the  cross  section  for  the  TRIGA  given  by 
West  (31)  the  value  of  the  disadvantage  factor  is  approximately 
1.37. 

(3.2)   Flux  Perturbations  Near  the  Control  Rods  and  the  Central 
Thimble 

For  proper  calculation  of  the  average  fltix  in  the  core  some 
correction  must  be  made  for  changes  in  the  neutron  flux  distribu- 
tion resulting  from  various  irregularities  in  the  core  lattice. 
During  normal  operation  of  the  Kansas  State  TRIGA,  one  control  rod, 
the  pulse-safety  rod,  is  completely  withdrawn  from  the  core,  but 
the  regulating  and  the  shim  rods  are  partially  inserted.   By 
assuming  that  there  is  a  flat  neutron  distribution  in  the  area  of 
these  control  rods  and  that  there  is  absorption  of  only  thermal 
neutrons,  flux  depressions  near  the  surface  of  the  rods  and  flux 
peaking  near  the  water  colximns  below  the  control  rods  and  in  the 
central  thimble  can  be  calculated  from  the  following  expression 
(14): 


where 


*  =  .^^fl  -  BK^(icr)/*  1  (21) 


<^  =  the  neutron  flux  in  the  vicinity  of  the 
control  rod, 

6     =  the  neutron  flux  in  an  unperturbed  region, 
o 

B  =  the  experimentally  determined  proportionality 
constant  for  the  effects  of  a  control  rod  on 
neutron  flux, 
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K   =  lAth' 

L   =  the  thermal  diffusion  length  in  the  core, 
th 

r  =  the  distance  from  the  center  of  the  control 

rod,  and 

K  (Kr)  =  the  modified  Bessel  function  of  the  second 
°      kind  of  zero  order. 

Experimental  data  obtained  with  the  miniature  fission  chamber 

can  now  be  used  to  calculate  values  for  B  at  various  vertical 

levels  near  a  control  rod.   Solving  equation  (21)  for  B  yields 

♦  -  ♦ 

By  using  flux  values  from  flux  insertion  hole  M  for  (t>  and  flux 
values  from  hole  H  for  *^,   experimental  B  values  are  obtained  for 
the  regulating  rod.   Although  similar  measurements  of  B  cannot  be 
done  for  the  shim  and  the  pulse-safety  rod,  it  is  possible  to  get 
a  good  approximation  for  B  for  these  rods,  assuming  (1)  propor- 
tionality between  B  and  the  rod  worth  and  (2)  water  columns  in  all 
three  rods  have  the  same  negative  B  values. 

(2.3)   Reactor  Power  Level 

An  IBM  computer  code  described  in  Appendix  (A-2)  is  used 
for  calculation  of  the  average  neutron  flux  in  the  fuel.   This 
program  uses  flux  data  obtained  with  the  miniature  fission  cham- 
ber and  accounts  for  both  the  disadvantage  factor  correction  and 
the  correction  for  flux  pertvirbations  near  the  control  rods  and 
the  central  thimble.   The  results  of  the  average  flux  calcula- 
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tions  and  the  values  for  reactor  power  level  calculated  from 
equation  (18)  are  summarized  in  Table  XI.   Detailed  numerical  cal- 
culations of  the  reactor  power  level  and  an  estimate  of  the  un- 
certainty associated  with  these  calculations  are  performed  in 
Appendix  {A-2), 
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(3.0)      Stnnmary  of  Results 

Figures  31  and  32  are  two  dimensional  contour  plots  summar- 
izing neutron  flux  data  taken  with  the  fission  chamber.   Figure 
31  shows  the  flux  contour  in  a  relatively  unperturbed  region,  and 
Fig.  32  shows  the  contour  on  a  plane  passing  near  the  regulating 
rod.   Table  XII  applies  to  both  Fig.  31  and  32.  The  contour  lines 
were  drawn  by  hand  with  no  attempt  to  interpolate  between  data 
points . 

These  plots  are  presented  in  this  manner  so  that  they  can  be 
compared  to  similar  plots  presented  by  West  (31).  Figure 
33  is  a  two  dimensional  plot  reproduced  from  West's  report  and 
applies  to  a  TRIGA  Mark  I  or  II  core  containing  56  fuel  elements 
and  four  completely  withdrawn  control  rods.   Table  XIII  applies 
to  this  figure.   This  plot  should  not  be  compared  directly  to  the 
experimental  data  for  a  number  of  reasons.   First,  the  neutron 
flux  in  Fig.  33  is  defined  by  (|)  =  nv  rather  than  by  "Westcott 
convention,"  (|>  =  nv^.   Second,  flux  in  Fig.  33  is  given  for  the 
homogenized  core  as  calculated  by  West;  whereas,  the  experimental 
data  are  taken  with  the  fission  chanber  in  the  water  gaps  between 
the  elements.  Third,  the  plot  from  West  is  for  thermal  flux 
(i.e.,  0  to  0.05  ev) ,  but  the  experimental  data  have  not  corrected 
for  epithermal  fission  of  U-235  in  the  chamber  and  does  not, 
therefore,  represent  strictly  thermal  flvix. 

Experimental  flux  values  measured  in  various  positions  in 
the  Kansas  State  TRIGA  are  summarized  in  Table  XIV.   The  values 
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reported  in  this  Table  are  calculated  by  using  the  beta-gararaa 
coincidence  method  for  calibration  of  absolute  flux. 
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Table  XII 

Neutron  Flux  for  the  Experimental  Two 

Dimensional  Contour  Plots 

(Figures  31  and  32) 

Westcott  Flux  Measured 

Contour 

Line              with  the  Fission  Chamber 

No. 

(n/cm^-sec-watt  x  10'') 

4 

0.500 

0.625 

6 

0.750 

7 

0.875 

8 

1.000 

9 

1.125 

10 

1.250 

11 

1.375 

12 

1.500 

13 

1.625 

14 

1.750 

15 

1.875 

16 

2.000 

17 

2.125 

18 

2.250 

19 

2.375 

20 

2.500 
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Figure  31.   Two  dimensional  flux  contour  of  the  TRIGA 
core  measured  with  the  miniature  fission 
chamber. 
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Figure  32.   Iwo  dimensional  flux  contour  of  the  TRIGA 
core  measured  with  the  miniature  fission 
chamber .  »•• 


70 


Table  XIII 


Neutron  Flux  for  the  Two  Dimensional 
Contour  Plot  Reproduced  From  West  (31) 
(Figure  33) 

Thermal  Neutron* 
Contour  Line  Flux 
No. (n/cm^- sec-watt  x  10'^) 

0  0.500 

1  0.550 

2  0.600 

3  0.650 

4  0.700 

5  0.750 

6  0.800 

7  0.850 

8  0.900 

9  0.950 

10  1.000 

11  1.050 

12  1.100 


*     Neutron  flux  as  used  by  West  is  defined  as  <|)  =  nv  rather 
than  by  "Westcott  convention'.' 
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Figure  33.   Two  dimensional  thermal  flux  contour 
as  presented  by  West  (31) . 
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(4.0)   Suggestions  for  Further  Study 

It  is  suggested  that  additional  work  be  performed  to  im- 
prove the  accuracy  of  the  absolute  flux  measurements  at  the  cali- 
bration point  in  the  TRIGA  thermal  colvimn.   Improvement  of  the 
beta-gamma  coincidence  technique  is  suggested  through  the  use  of 
a  high  efficiency  beta  detector  and  more  sophisticated  pulse 
shaping  circuitry.   Further  development  of  the  "s\im  peak"  method 
is  also  suggested,  with  particular  emphasis  on  methods  for  cal- 
culating the  areas  under  the  photopeaks . 

For  this  work  thermal,  resonance,  and  fast  flux  measurements 
were  performed  only  in  the  central  thimble  and  the  rotary  speci- 
men rack.   Similar  measurements  could  be  made  for  the  other  irra- 
diation facilities  of  the  TRIGA  Mark  II  (i.e.,  the  pneumatic 
transfer  system  and  the  four  beam  ports) . 

The  detailed  experimental  flux  plots  obtained  using  the 
Reuter-Stokes  fission  chamber  may  be  compared  with  the  neutron 
flux  calculated  by  various  analytical  methods.  When  the  IBM-360 
computer  is  installed  at  Kansas  State  it  should  be  feasible  to 
perform  large  multi-group  calculations  to  obtain  flux  plots 
similar  to  those  presented  by  West  (31) .   It  might  also  be  feas- 
ible to  devise  a  two  dimensional  code  which  calculates  the  flux 
perturbations  near  the  control  rods  and  the  other  irregularities 
in  the  core.   These  computer  codes  could  then  be  used  to  calcu- 
late the  flux  shapes  for  the  actual  core  loading  in  the  Kansas 
State  TRIGA  for  various  control  rod  positions  and  compared 
directly  to  the  experimental  flux. 
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APPENDIX 

(A-1)   Weighted  Least  Squares  Analysis  for  Calculating  Initial 
Count  Rate 

The  count  rate  at  time  t  of  a  decaying  radioisotope  is 
given  by  the  relation 

InC  =  -Xt  +  InC  (A-1) 

where 

C  =  the  count  rate  at  time  t, 

C^  =  the  initial  count  rate, 

X   =  the  radioactive  decay  constant  of  the 
isotope,  and 

t  =  the  time  after  the  initial  count  rate. 

The  program  described  here  fits  the  linear  relation  of  equation 

(A-1)  to  a  series  of  count  rates  using  a  weighted  least  squares 

method.   The  error  associated  with  each  count  rate  is  given  by 

*i  "   5-?f^"^i  "  ^-^^i  +  lnC^)l-      (A-2) 
In  this  method  the  sum  of  the  errors  squared  is  written  as 

n      n 

^'Tj'I'  Z)^?^^^^i  ^  '^i  -  ^"^o>'  <^-3) 

i=l     i=l  ^ 
where 

n  =  the  number  of  count  rates  in  the  fit, 

1  th 

-2  =  a  weighting  factor  for  the  i   count  rate, 

i 
a^   =  the  variance  on  the  count  rate  after  cor- 
rection for  background  and  dead  time. 

Since  A  is  a  known  value,  it  is  held  constant  and  the  least 

squares  criterion  is  satisfied  by  choosing  InC  such  that  E  is 
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a  minimum.   This  minimization  is  done  by  taking  the  derivation 

of  E  with  respect  to  InC   and  setting  this  derivation  equal  to 

zero  as  follows: 

n 
dE 


ddnC^)  -  °  - 


=  2y^^2(lnC.  +  At.  -  InC  )  (-1)  {A-4) 

/  J  o .      1       1        o 


i=l 


Solving  for  the  InC  yields 

n 


Z-2  (At.  +  Inc.  ) 
of    1       1 


mC^  =  i^^^ ,  (A-5) 


T^i 


1=1  ■ 

and  the  initial  count  rate  is  found  by  simply  taking  the  anti- 

logrithm  of  equation  (A-5) .   The  variance  of  the  initial  count 

rate  C   is  calculated  from  the  expression 

n 


C  2y^i2(lnC   -  At.  -  Inc.) 

o  /  yg. '     O       1        1 


ol     =    ^'^  \ .  (A-6) 

o 


[n-e 


2 

i 
i=l 

As  a  check  to  determine  whether  or  not  the  isotope  is  de- 
caying at  the  expected  half  life,  the  program  also  performs  a  two 
parameter  fit  using  the  least  squares  technique  and  solves  for 
the  radioactive  decay  constant  A .   In  addition  to  the  least 
squares  analysis,  the  program  corrects  for  the  background  count- 
ing rate,  the  dead  time  of  the  counting  system,  and  the  decay  of 
the  isotope  during  a  finite  counting  time.   The  dead  time  calcu- 
lated for  the  Geiger-Mueller  counting  system  described  in  section 
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(1.2)  was  172.8  micro  seconds.   The  final  result  of  the  program 
is  the  initial  saturated  count  rate  in  units  of  counts/minute- 
gram. 
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Fortran  Program  for  Calculation  of  Initial  Counting  Rates 


^'CNi.tS      JOB   INITIAL  ACTIVITIFS 

MCN5$      CCf-'T  15  N^INUTES?  1  .■  PAGF  •    UCHEY*  N.  -.  ^rnr. 

MCNIS      ASGN  MJB»12 

MCN'f.$      A5GM  MGC»16 

VCNIX       iCDE  GO  »  TEST 

'■■^CNiS      EXEQ  FORTRAN,  ,,»,  ,  .EXTRAPU 

0IK.ENS1CNAC(8)  »  VvF  (  10  )  ,  T  I  LPL  I(  10  )  ,  l)  (  3  ) 

1  FORMAT{2E20.8) 

2  FORMAT  (  I1^,3F1C.1,  11..) 
I:    '^0RMAT(8A1C) 

1   ^0RMAT(6F1C. 1) 

]]  F0PMAT(24H!<INDICATrD  HALF-LIFE  IS»E14.8»5H  MIN.) 

]?  FORMAT  (29HT  I  NI  T  ^ '^i   cnr^ycir  ^CUNT  RATf"=  fC-'' '^  .° -^u  rp'  /r/- .  ^  X  ,  ! '^'HDEV  I 

?ATION'=,FTA.a) 
1^  F0RM„T(60HTINITI AL  COL'  '    TtS  NORMALIZED  TOINFINITE  IRRADIATION 
2TIKE) 

READ( 1 ,1 ; :  -  . ,TAU2 

READ( 1 »2) NFOIL 

-VRITE(3»13) 

D03H=1 ,NFOIL 

READ(1,4) AC 

WRITF(3,4)AC 

PFADd  ,2)  N'TC,TIR,FVASS,TC»MAT 

SEC=0. 

5FCT  =  :". 

?T  =  i.. 

ST?  =  i  . 

SWF  =  i.'. 

SESG=C. 

SFNLG=^ . 

GOTO  (3  »6.7)  f.v.AT 

5  A=AL0G(2.  )/(  2.78*24.^*6;}.  ) 
G0T08 

6  A  =  AL0G(2.  )/54.:' 
GOT  08 

7  A=AL0G(2.  )/( ]4.3*24.*6r. ) 

8  TA  =  ALOG{  A*TC/(  1  , -EXP  ( -A*TC  )  )  )  / '^ 
PFAD(  1  »10  )  n(  1  )  ."^t  ?)  ,n(  ^  ) 
D09K=1 »NTC 

RFAD(  1  »1.'  )CO,Tn,Tn>TM,B0,TB 

T=6o,*(24.*TD+TH)+TM+TA 

C  =  C.  /(  TC-C^  *TAU)-r-,:/  (TB-BO*TAU) 

Ct.F  =  ]  .  /  (  TC-CC*TAU)**2 

[•;!.F  =  1./  (  TcJ-B0*TAU)**2 

'■  F(K)  =  1./D(i<)**2 
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]4 


EC=ALCG(C) 

SFC=?FC+FC 
SFCT  =  SECT-^P^*T 

5T=ST+T 

SWF  =  SWF  +  aF(,<) 

TILPLI (K) =A*T+FC 

SFNLC=SFNLO+WF(K)*TILPLI (K) 

RMTC=NTC 

PL(-  =  '^FNLO/SWF 

A(.  =  EXP(RL>  )/(FMASS*(  1  .-EXP  ( -A*T  I  R  )  )  ) 

THALF  =  -ALCG{2.  )  *(  RNTC*'' ^'^- "  T^  5T  )  /  (  RMTC*. 

'•.RITF(3,11  )THALF 

'^•Cl^f<  =  l»MTC 

5FS0  =  StSQ+'A'F  (.<)*(kLo-TILHLI  (.'.  )  )^-^2 

TDFV  =  A^*SQRT  (SESQ/  (  Sv-F*  (  RNTC-l  .  )  )  ) 

WRITE(3,12)A(J,TDEV 

STCP 

MCNSS  EXEQ    LINKLCAD 

CALL    EXTRAPd 

•••^OM  i  FXFO    FXTRAPB*'  .  _ 


rCT-'  T' 


T) 


?  -  E  E  ♦:•  •  ^""""'^ 
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(A- 2)   Calculation  of  Reactor  Power  Level  by  the  Integrated  Flux 
Method 

To  calculate  the  power  generated  in  the  TRIGA  core  the 
average  flux  in  the  fuel  portion  of  the  core,  ♦,  is  determined 
from  the  thermal  flux  data  taken  with  the  miniature  fission 
chamber.   To  calculate  this  average  flux  the  core  is  divided  into 
annular  elements  as  shown  previously  in  Fig.  30.  The  average 

flux  is  then  approximated  by  the  expression 

14   N 

i-1  j-1 
where  i  indicates  the  vertical  level,  ♦^j,  is  the  average  flux  in 
the  ij^^  volume  element,  and  H  is  the  number  of  fuel  elements. 
The  average  flux  in  a  volume  element,  ♦^j,  is  obtained  by  estimat- 
ing from  the  geometry  of  the  core  and  from  the  measured  values, 
the  flvix  at  each  of  the  four  corners  of  the  annular  element  and 
by  taking  the  average  of  these  values  divided  by  the  disadvantage 
factor.   The  computer  program  presented  at  the  end  of  this  appen- 
dix corrects  the  flux  data  for  flux  perturbations  near  the  cen- 
tral thimble  and  the  three  control  rods,  and  performs  the  numeri- 
cal integration  described  above.  Flux  data  taken  at  eight  dif- 
ferent radial  positions  (foil  insertion  holes  B,  C,  F,  G,  H,  K, 
O,  and  the  central  thimble)  are  used  in  the  program.   These  posi- 
tions were  chosen  because  they  represent  positions  where  control 
rod  effects  are  relatively  small. 

The  values  of  B  used  for  the  flux  perturbation  calculations 
are  given  in  Table  XII.   The  B  values  for  the  regulating  rod 
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were  determined  experimentally  by  using  equation  (22)  with  data 

from  foil  insertion  hole  H  as  unperturbed  flux  values  and  data 

from  foil  insertion  hole  M  as  perturbed  flux  values.  The  B 

values  for  the  shim  and  the  pulse-safety  control  rods  and  for 

the  central  thimble  were  estimated  from  the  B  values  calculated 

for  the  regulating  rod. 

Table  XII 

saaa>»snuaBsss:^aaBEccsaE<BS9sa«saasxsi3iss39Basxrsssast3)ss8aBiiss3cassisB:x3rssessssa>x:a;Bcs3ss>iss3^^atsi 

Calculation  of  the  Experimental 
Constant  B  in  the  Equation 
i   -  *^(1  -  BK^Cicr)/*^) 

Vertical  Regulating  Rod  Shim  Rod  Pulse-Safety  Rod  Central 
Level                                          Thimble 
(X  10^)        (X  loT) (X  10^) (x  10^) 


1 

0.77 

2,00 

-0.77 

-0.77 

2 

0,95 

2,25 

-1.35 

-1.35 

3 

0.95 

2,45 

-1,71 

-1.71 

4 

1,15 

-0.55 

-1,90 

-1.90 

5 

1,15 

-1,75 

-1.90 

-1.90 

6 

0,59 

-1,70 

-1.90 

-1.90 

7 

-0,55 

-1.90 

-1.90 

-1.90 

8 

-1,53 

-1.90 

-1.90 

-1.90 

9 

-1,54 

-1.90 

-1.90 

-1.90 

10 

-1,92 

-1.90 

-1.90 

-1.90 

11 

-1,90 

-1.90 

-1.90 

-1.90 

12 

-1.71 

-1.71 

-1.71 

-1.71 

13 

-1,35 

-1.35 

-1.35 

-1.35 

14 

-0.77 

-0.77 

-0.77 

-0.77 

As  shown  by  equation  (18)  the  power  generated  in  a  reactor 
is  simply  the  fission  rate  in  the  core  times  the  volume  of  fuel 
in  the  core  times  the  energy  released  per  fission.   However,  for 
the  calculation  of  power  it  is  convenient  to  rewrite  equation  (18) 

as 

KN  mo- 

A235 
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where 


m  =  the  mass  of  the  U-235  in  the  core, 
A235  =  the  atomic  weight  of  U-235, 

0^=  the  microscopic  fission  cross  section, 
N  =  Avogadro's  number,  and 
K  =  the  energy  released  per  fission. 
Substituting  the  proper  nxomerical  values  into  equation  (A-8) 

yields: 

(2.88  X  10"^M  (0.6023  x  10^'*)  (2325)  (580  x  lO"^'*)  (1.03  x  10^^ 
P  =  (235) 

P  =  1.03  watts, 
where  the  value  for  ?  was  obtained  using  the  beta  gamma  coinci- 
dence method  for  determining  absolute  flux  and  has  been  corrected 
for  the  ratio  of  the  maximum  flux  in  the  water  to  the  average 
flux  in  the  fuel.  The  microscopic  fission  cross  section  for 
U-235  at  2200  meters  per  second  was  tsOcen  as  580  *  7  bams. 
Factors  for  correcting  the  cross  section  for  departure  from  the 
1/v  form  are  tabulated  by  C.  H.  Westcott  in  Effective  Neutron 
Cross  Sections.   Since  this  pxiblication  was  not  available  the 
2200  meters  per  second  cross  section  was  assumed.   It  is  estimated 
that  this  correction  would  lower  the  effective  cross  section  by 
somewhat  less  than  2%.   It  is  not  necessary  to  correct  the  cross 
section  for  the  epithermal  component  of  the  neutron  spectrum 
because  the  flux  data  was  not  corrected  for  epithermal  fission 
in  the  fission  chamber. 

To  obtain  an   estimate  of  the  inaccuracies  in  the  power 
calibration  each  of  the  quantities  in  equation  (A-8)  are  consid- 
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ered  separately.   The  values  for  Avogadro's  nvunber,  the  mass  of 
the  fuel  present  in  the  core,  and  the  atomic  weight  have  neglig- 
ible uncertainty  compared  to  the  other  quantities.   The  total 
nuclear  energy  released  per  fission  is  (3.16  +  0.13)  x  10  * ^ 
watt-sec. /fission.   Additional  uncertainty  results  in  this  con- 
stant when  estimating  the  actual  calorimetric  power  output.   The 
value  of  K  in  this  case  becomes  (2.88  +  0.20)  x  lO"^! 
watt-sec. /fission.   The  microscopic  fission  cross  section  as 
already  stated  is  580  ±   7  barns.   The  final  uncertainty  in  the 
power  calculation  results  from  the  difficulty  in  calculating  the 
average  fl\ix  in  the  fuel.   An  estimate  of  the  statistical  devia- 
tion on  the  average  flux  in  the  fuel  can  be  obtained  from  the 
following  expression: 


m  M 


o2  =  Ly^  y^  v^a?  .  (A-9) 


where 


ik 
i=l  k=l 

a^  =  the  variance  on  the  average  flux  in  the 
"i       fuel, 

(fr.,  =  the  flux  value  measured  with  the  fission 
chamber , 

a.        =  the  variance  on  <{>.,,, 
*ik  ^^ 

V  =  the  total  volxime  of  the  fuel, 

V,  =  the  volume  of  fuel  effected  by  (ji., 

M  =  the  number  of  flux  measurements  per  verti- 
cal level. 

Only  eight  independent  flux  measurements  were  used  across 
the  core  at  each  vertical  level  in  the  average  flux  calculation 
euid  all  other  values  of  <fr.,  were  estimated  from  these  eight 
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measurements.   (See  page  83)   It  follows  then  that  in  calculating 
the  statistical  variance  on  the  average  flux  from  equation  (A-9) 
that  the  summation  should  be  taken  only  over  the  number  of  inde- 
pendent measurements  (i.e.,  k  runs  from  k  =  1  to  k  =  8  at  each 
vertical  level) .   V,  is  a  factor  which  weights  the  flux  by  the 
volume  of  the  fuel  estimated  to  produce  that  flux.   Vj^  is  esti- 
mated as  shown  in  the  simplified  sample  calculation  given  below. 

4)12  ♦is 


In  this  example  (\>i,    (jj^,  and  (pn    indicate  fliix  values 
measured  by  the  fission  chamber  with  associated  variances  given 
by  a .  ,  a ,  ,  and  a ,   .   The  average  flux  in  the  fuel  elements  of 
the  sketch  is  calculated  as  discussed  previously  by  assuming  that 

(j)^  =  (|,2  =  (j)3, 

and  ^ 

Equation  (A-9)  for  calculating  the  variance  on  ^  in  this  example 
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sums  only  over  the  three  measured  flux  values  (^i,  ♦!,/  and  ^u). 
The  volume  weighting  factor  V.  for  each  of  these  ineasureiaents  is 
determined  by  estimating  a  volume  of  fuel  associated  with  the 
points  at  which  the  flux  is  equated  to  tlie  measured  flux.  For 
example  to  deto3rmine  Vj ,  the  weighting  factor  for  <^x,  it  is 
assumed  that  h   of  fuel  element  number  1  is  associated  with  ^i, 
that  H   of  fuel  element  1  and  H   of  fuel  element  2  is  associated 
with  ^2 1   an^  that  h   of  fuel  element  2  is  associated  with  ♦g. 
Thus,  Vj ,  represented  by  the  louer   cross  hatched  area  in  the 
sketch,  is  equal  to  the  arm   of  these  volumes  or  1  volume  unit, 
where  each  circular  fuel  element  represents  a  xinit  voluine. 
Similarly  V^,  the  weighting  factor  for  <>i»  represented  by  the 
\incross  hatched  fuel  elements  in  the  center  of  the  sketch,  is 
equal  to  2.5  volxmo  units.  Vji,  represented  by  the  upper  cross 
hatched  elements  is  equal  to  1.5  volume  units,   now  suljstituting 
into  equation  (A-9)  the  following  result  is  obtained 

o2  s  i^ial      +  6.25  al     +  2.25  a^      ).       (A-10) 
•^        25   ♦x  ♦!»  ♦!! 

The  variance  calculated  from  this  equation  is  ±  5.6%  of  ^. 
Additional  error  in  calculating  the  average  flux  in  the  fuel 
results  from  the  uncertainty  in  calculating  the  perturbations 
near  the  control  rods  and  the  central  thimble.  Extensive  calcu- 
lations of  this  sort  would  require  a  large  computer  code  and 
would  probably  not  be  practical  on  the  IBM  1410  at  Kansas  State; 
therefore  the  approximation  of  equation  (21)  was  used.  Errors 
in  calculating  the  flux  perturbations  from  this  expression  may 
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be  as  large  as  ±  50%,  but  the  effect  of  the  local  perturbations 
changes  the  total  power  level  by  an  estimated  +  5%,  making  the 
uncertainty  on  the  total  power  level  due  to  local  perturbation 
±  2.5%.   Assuming  that  errors  due  to  the  correction  for  local 
perturbations  are  independent  of  the  statistical  errors  and  com- 
bining these  tv/o  values,  the  value  of  average  flux  in  the  fuel 
becomes  (1.03  ±  0.07)  x  10^  n/cra^-sec-watt.   Finally,  after 
combining  all  the  uncertainties  the  calculated  power  level  for  a 
power  of  one  watt  indicated  by  the  linear  recorder  of  the  reactor 
console  is  1,03  ±  0.11  watts. 
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FORTRAN  Program  for  Calculation  of  the  Average  Flux  in  the  Fuel 


O  w  i->     '  v^  v:  ..  1  -^   \_  '-'  i_  i  -   / A  I  i  w  I  . 

r'CN-S      CCMT  15  MINUTES*  10  PAGf   •    UCHEY  ^   .   .  OtPT. 

MCNSS      CG.MT  USPS  PHYSICAL  TAPE  '^ 

N'ON?-^      ASGN  MJB»12 

MCr-'T-J      ASGN  MGC»16 

MCN$$      MODE  GCjTEST 

MCNIS      EXEO  FGRTRAN»» » .» » »PCALId 

DIiVENSICNAC  10,30)  ,RSH(8)  ,RR  (  8  )  » RP  (  8  )  ,RCT  (  8  )  ,SUM(5)— 

CCM.'-ICNRLTH 

1  FORMAT {8E1C, 3) 

2  FCRMAT(5X,13HAVERAGE  FLUX= IPE 14  .  7 ) 
?  FORMAT (8F10. 1 ) 

PLTH=1.6  ' "^" 
FLLIX='':. 

Z  =  :  . 

PFADd  ,3)RSH 
"FADd  ♦3)RR 
READ{1,3)RP 
READd  ,3)RCT(  1  ) 
D016I=2.6 
16    RCT( I )=RCT{ 1 ) 
D027J=1,14 

R  EAD ( 1 , 1 ) A ( 1 , 1  )  , A ( 2  » 1)  » A ( 3  >  1  )  ,  A ( 4  ,  i )  , A ( 5 , 1 )  » A ( 6  » 1 )  >  A ( 8  » 1 )  » A (  1  0 , 2 1 ) 
D06I=2»6 

6  A  (  1,I)=A( 1,1  ) 
D07I=?,6 

7  A(2,I)=A(2,1) 

nrei=?,i2 

8  A(3,I )=A( 3,1) 
D09I=2,12 

9  A(4,I )=A(4,1 ) 
DOK  1=2,18 

]       A  (^,1  )=A(  5,1  ) 

D011I=2,18 
11    A(6,I )=A(6,1 ) 

00121=2,24 
1?    A  (7,1  )=A(6,1  ) 

'^r]3I=?,24 
I''     A(8,I)=A(8,1) 

D014I=2,3w 
Ih    A(9,I)=A(8,1) 

D015I=2?,27 
15    Ad  :  ,1  )=A(  10,21  ) 

READd  ,1  )  BR^-'^M.RP 

BCT=RP 
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RFWIND5 

V.'RITE(5)  (  A(  1,1  )  ,1=1  ,6)  .7.7 

CALLCCRR(RCT,BCT»6) 

RFAD(5 )  ( A( 1,1 )  ,  1  =  1  ,6  ) 

REWIND3 

WRITE(5)A(2,2)  ,A(3,3)>A(h,2)»A(4,3)»A(5,3),A(3,4)  ,Z,Z 

CALLCCHR(RSH,bSH,6) 

READ(5)A(3,2),A(3,3)»A(4,2),A(4,3)»A(5,3),A(5,4) 

REWIND5 

.•.'RITE(5)A(4,6)  ,A(4,7  )  ,A(5,9)  ,A(  5,  K  )  ,A(6,9)  »A(6»10)  ,A(7,12)  ,A(7,13 
B) 

CALLCCRR( RP,BP»8) 

RFAD(5)A(4,6),A(4,7),A{5,9),A(5,1C),A(6,9),A(6,10),A(7,]2),A(7,13) 

RFWIND5 

WRITE(f3)A(6,15),A(5,16),A(7»20)»A(7,21),A(8»20)»A(8,21),A{9,26),A( 
B9,27) 

CALLCQkR{Kk,D,<,8) 

READ(5)A(6»15)»A(6,16),A(7,20)»A(7»21),A(8,20),A(8,21),A(9,26)>A(9 
B,27) 

L  =  l 

V  =  l 

DC25<=5,24,6 

SUN'  (  L  )  =  (  A  {  M  ,  1  )  +  A  ( '•'+ 1  ,  <)  +  A  (  V  ,  K  )  +  A  (  M+  1  ,  1  )  )  /4  , 

DC241=2,K 
2  4  SUM  (I  )=SUM(L)  +  (  A(^',  I -1  ) +A  (  M  ,  I  ) +A  (  M  +  1  ,  I-l  ) +A  (  M+ ]  ,  I  )  )  /4. 

L  =  L  +  ] 

SUM ( 2 ) =SUM ( 2 ) - ( A ( 3 , 2 ) +A ( 3 , 3 ) +A ( 4 , 2 ) +A ( 4 , 3 ) ) /4 . 
SUM  {  3  )  =SUM  (  3  )  -  (  A  (  b  ,  9  )  +A  (  5  ,  1 0  )  +A  (  6  ,  "^  )  +A  (  6  ,  10  )  )  /4. 
SUM ( 4 ) =SUM ( 4 ) - ( A ( 7 , 2b ) +A ( 7 , 2 1  ) +A ( b  ,  20 ) +A { 6 , 2 1 )  ) /4 . 
SUM(5)=-( A(9,24)+A(9,2  5)+A( IC ,2  4)+M( 10,2  5) )/4. 
DC26I=22,26 
2  6  SUM(  5)=SUVi(5  )  +  (  A(9,  I  )+A{Q,  I  +  l  ) -t- a  (  IC,  I  ) +A  (  10  ,  I  +  1  )  )/4. 

DC27N=1 ,5 
27  FLUX=FLUX+SUM(N) 
VRITF(3,2)FLUX 
STCP 
END 
MCNSS      EXEO  FORTRAN 
SUBRCUTINECCRR(R»B,N) 
DIMENSICNP(8) ,R(8 ) 
CCMMCNRLTH 
REWINDS 
READ(5)P 

'?EWI^'D5  .!>  .. 

PCI  1  =  1  »N        -^  ••   \ 

nFS  =  RESK(C,R(  I  ) /RLTH)  "7^, 

IFOFS.LT.C.  )«ES  =  0. 
1  P (  I  )=P( I) *( I .-B*PFS) /P (  I ) 
»'RITF(5)P 
REWIND5 

RETURN 
bND 
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••V?^"•t  3  EXEQ    FORTRAN 

'^Ll^'CTIC^'    APPSER(X»N»SIGN'/i  »  T  ) 
IF( I.EO.O)GC    TC    2 

TEPM=(4.*FLCAT(N**2)-1. )/( B.^X) 
ANPLUS=3. 
AN  =  2. 
GO    TC    3 
TERM=1. 
ANPLUS 
A\'=l. 
ApPSFR 
^vijLT  = 
5IGNA= 

X  =  APP<^ 

TERKE= 
.    TERf^  =  T 
1.  )  ) 


=  1  . 

=  JFPV1 

( (8.*X)**2)*?IGNA 

4.«-FLCAT(N**2) 

FR 

TERM 

ERM*(SIGNA-ANPLUS**2  )*(  SICiNA- (  ANPLUS+2  .- )  **2  )  /  (  A--'ULT*AN*  (  AN+ 1 


21 
22 


12 

15 

3 


IF(APS(Tlk'-;c  )  .  Ll./<l'S  (  TERr  )  )' 

APPSER=APPSER+TERM 

IF(X.Ew.APPSER)GC  TC  7 

X=APPStR 

TERME=TERM 

AN=AN+2. 

ANPLi'<:  =  ANPLU5+4. 

GC    TC  A 

RETURN 

'CNri  EXEQ    FORTRAN 

FUNCTION    dESSER(N»X»SI6NA, I ) 
GAN'MA=.5772156649C 
IF(N.Nt.u )G0    TO    2 
SUf'  =  l. 
GO    TC    22 
SU''1  =  1./FACT(M) 
DC    21    J=l ,N 
5UN'  =  SUM*(  X/2.  ) 
A  ^lPLl  |C;  =  \i+l 

I F(  I  .PQ.r  )GC    TC    3 

VFRM^ALCG(X/2. )+GAmmA 

IF(N.EQ.0)GO    TC    15 

DC    12    J=1,N 

VERK  =  VERM-.5/FL0AT(  J ) 

SU>  =  SUf''i*VERM 

TERiM  =  SUK 

Fr/ULT  =  X*X/4.*SIGNA 

X  =  l. 

sur-*E=.suv 

IF( I.FG.O )G0    TC    11 
TER'>^  =  TERM/VERf/ 
VERM=VERM-.5/X-.5/A\'PLU.S 
TERM  =  TERM*VERN' 


TC    7 
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]^     TFP-X  =  TtRM*FMULT/(ANPLU5*X) 

5U^»  =  SUV  +  TER^^ 

6  ANPLllS  =  A'\|PLU<^+i  . 

X=X+1. 

IF{SUME.NE.SUi-')GC  TC  4 

7  BESSEP  =  c'J-'^- 

RETUR^, 

LNL- 

v:OMj>S>      iiXt-U  Fi.kJ :\Ai\ 

FUNCTION  BESSES(N»X,SIGNA) 

RFSSFS=0. 

IF(N.LF.C)GC  TT  1 

no  7  J=1,V 

vfM  =  -Ki+2*{  J-1  ) 

IF(MM.EQ.':  )GC  TC  A 

!F(MM.GT.C)GC  TC  5 

AMULT=(2./X)*-*(-MN) 

G  C  1  \^    6 

4  AMULT=1. 

GC  TC  6 

5    AMULT= (X/2. ) **MN 

6  BESSES  =  BESStS+(  (  S  IGNA  )  **  (  J  + 1  )  )  *FACT  ( ,\-J  )  ^A.-iULT/F  ACT  (  J-1  ) 

7  CCNTINUE 

1  RFTURK 

FMD 

^r^i':^      FXFO  FORTRAN 

FDNCTICN  BFSK{N,X) 

1   FCRVATdH  ,27HCRDER  N  CF  BESK  IS  NEGATIVE) 

1;  PCRMATCIH  ♦3CHARGUMENT  X  CF  BESK  IS  NEGATIVE) 

12  FCRMATdH  ,33hARGUMENT  X  CF  3ESK  lo  TCC  LARGE  »  =  ♦£  14 .  o  ) 

PI  =  3.1415  926  53  5  8979  3?''84626 

IF(N.G£.0 )GC  TC' 2 

r:=-N 

v.RITE(3.K) 

2  IF(X.EQ.O.)GC  TC  9 

IF{X.GT.O.)GC  TC  4 

X=-X 

VPITF(3,11) 

4  IF(X.LT.227.5)GC  TC  6 

BESK=0. 

WRITEC3»12)X 

RETURN 

6  IF{X.G^.^.3)GC  TC  6 

bESK  =  .5*B£SSES(>J»X,-l.  )  +  ((-I.)«*(i^+l)  )  ^liESbER  (  N  ,  X  ,  1  .  ,  1  ) 

RETURN 

P  BESK  =  SC.RT(PI/(2.*X)  )*EXP(-X)«  (  APPSER  (  X  ,  N  ,  1  .  »  0  ) +APP5ER  (  X  ,N  ,  1  . 

»1)  ) 

RETURN 

9  BFSK  =  .999999999999999999E  +  9'^; 

RFTURN 

FND  . 
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7  F 
F 
I 

1  F 
I 

3  F 
R 

4  \j 
F 
C 
I 

5  C 

2  R 

6  W 


UNCT 
ORMA 

ACT=: 

FIN) 
ACT  = 
F(N- 
ACT  = 
ETUR 
C  5 
ACT  = 
ALLO 
F(  J. 
CNTI 
FTUR 
RITF 
ETUR 

:nd 

;f\!i;  i 
;  N 1  i 


IC 
T( 
1. 

1» 


EXEC  FORTRAN 
N  FAC.T(N) 

IH  ♦3CHARGUMENT  F; 


■R  FACT  II,  T; 


LARGF) 


2»4 

9999999  999999.9999E9  9 

/2)*2)2»3»2 

ACT 


(\ 

-F 

N 

1  = 

FA 

VF 

EQ.l )GC  T 

NUE 

N 

(3.7) 


1,N 

CT*FLOAT( I ) 
RFL( J) 


EXEG  LINKLCAD 
CALL  PCALIB 

EXEG  PCALIS,MJb- 
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NOMENCLATURE 

A 

activity,  (dis./sec.) 

^. 

activity  of  a  bare  foil 

^cd 

activity  of  a  cadmium  covered  foil 

\ 

atomic  weight  of  phosphorous 

^s 

atomic  weight  of  sulfur 

^ 

zero  time  activity 

^ 

area  under  photopeak  #1 

^2 

area  under  photopeak  #2 

^12 

area  under  the  sum  peak 

^235 

atomic  weight  of  U-235 

B 

experimentally  determined  proportionality  constant 
for  the  effects  of  a  control  rod  on  the  neutron  flux 

C 

count  rate 

% 

initial  saturated  count  rate  from  the  phosphorous 
pellet 

^8 

initial  saturated  count  rate  from  the  sulfur  pellet 

^O 

zero  time  count  rate 

E 

energy 

«1  ^  ^2 

probabilities  that  gamma  1  and  gamma  2  are  absorbed 
in  the  crystal  by  a  photoelectric  process 

^cd 

correction  for  epithermal  absorption  of  neutrons 
in  cadmivmi 

^h 

correction  for  flux  depression  in  the  vincinity  of 
a  foil 

K 

energy  released  per  fission 

^8 

32 
isotopic  content  of  S 

K^(<r) 

zero  order  modified  Bessel  function  of  the  second 
kind 
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m 

mass 

"a 

Avogadro's  number 

n 

total  neutron  density,  (n/cm  ) 

% 

the  background  coincidence  count  rate 

""r 

random  coincidence  count  rate 

"^e 

count  rate  from  the  beta  detector 

'^Y 

count  rate  from  the  gamma  detector 

"eY 

beta-gamma  coincidence  count  rate 

-h 

total  observed  beta-gamma  coincidence  count  rate 

p 

power 

R 

cadmium  ratio 

r 

radial  distance 

"p 

content  of  P^^  in  the  phosphorous  pellet  in  weight 
percent 

T 

total  area  under  the  gamma  spectrum 

t 

time 

H   ^  ^2 

probabilities  that  gamma  1  and  gamma  2  are  absorbed 
in  the  crystal  by  any  means 

V 

volume 

^c 

volume  of  the  fuel  in  the  core 

V 

mean  velocity  of  neutrons  in  the  reactor  core 

^o 

most  probable  neutron  speed  in  a  Maxwellian  distri- 
bution at  293°  K 

a 

ratio  of  the  resonance  component  of  the  resonance 
activation  to  the  1/v  component 

h 

4"H 

least  squares  error  associated  with  the  i   coiint 
rate 

^3^S 

efficiencies  of  the  beta  and  gamma  detectors  res- 
pectively 
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K         ^Ath'  w^«^e  ^th  ^^  ^^®  thermal  diffusion  length 

X         radioactive  decay  constant 

Z         macroscopic  absorption  cross  section 


a 


I         macroscopic  fission  cross  section 

0^         standard  deviation  on  activity  A 
A 


'=o 


standard  deviation  on  the  coiint  rate  C 


o(E)       microscopic  activation  cross  section  as  a  function 
of  energy 


0         microscopic  fission  cross  section 

f 

•i 
'p 


th 
standard  deviation  on  the  i   count  rate 

0         microscopic  cross  section  for  P 


its   dE/E   resonance  portion  of  the  total  resonance  activation 
•^  ^  integral 

la,.   dE/E  1/v  portion  of  the  total  resonance  activation  integral 

'^  32      32 

saturation  cross  section  for  the  S   (n,p)P   reaction 
o 

standard  deviation  on  the  flux  ^ 

resolving  time  of  the  coincidence  circuit 

2 
neutron  flux,  (n/cm  -sec) 

integrated  fast  flux  above  2.9  Mev 

resonance  neutron  flux 

r 

thermal  neutron  flux 
th 
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The  purpose  of  this  work  is  to  determine  experimentally  the 
neutron  flux  in  the  Kansas  State  University  TRIGA  Mark  II  nuclear 
reactor  and  to  use  the  measured  flux  to  calculate  the  power  gen- 
erated by  the  reactor.   "Westcott  convention"  was  used  throughout 
this  work  for  reporting  neutron  flux  because  of  the  convenience  of 
this  definition  for  experiments  involving  the  calculation  of  re- 
action rates.   Absolute  thermal  neutron  flux  is  measured  at  a 
selected  position  in  the  thermal  column  and  used  as  a  standard 
for  calibration  of  all  other  flux  measurements.   The  calibration 
point  in  the  thermal  column  also  provides  a  known  thermal  flux 
which  may  be  useful  for  many  experiments  at  Kansas  State.  Ther- 
mal, resonance,  and  fast  fluxes  are  measured  at  the  point  of 
maximum  flux  in  the  central  thimble  and  in  the  rotary  specimen 
rack.  Vertical  plots  of  thermal  flux  are  obtained  with  a  minia- 
ture fission  chamber  inserted  into  the  reactor  core  through  the 
foil  insertion  holes  in  the  top  grid  plate.   By  numerical  inte- 
gration of  the  thermal  flux  data  obtained  with  the  fission  cham- 
ber, the  average  flux  and  the  corresponding  power  output  of  the 
core  are  calculated.   At  an  indicated  power  level  of  one  watt  on 
the  reactor  console,  the  power  calculated  from  the  flux  data  was 
1.03  ±  0.11  watts. 


